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Summary 

During the SARS-CoV-2 pandemic, novel and traditional vaccine strategies have been deployed 
globally. We investigated whether antibodies stimulated by mRNA vaccination (BNT162b2), 
including 3™ dose boosting, differ from those generated by infection or adenoviral (ChAdOx1-S 
and Gam-COVID- Vac) or inactivated viral (BBIBP-CorV) vaccines. We analyzed human lymph 
nodes after infection or mRNA vaccination for correlates of serological differences. Antibody 
breadth against viral variants is less after infection compared to all vaccines evaluated, but 
improves over several months. Viral variant infection elicits variant-specific antibodies, but prior 
mRNA vaccination imprints serological responses toward Wuhan-Hu-1 rather than variant 
antigens. In contrast to disrupted germinal centers (GCs) in lymph nodes during infection, 
mRNA vaccination stimulates robust GCs containing vaccine mRNA and spike antigen up to 8 
weeks post-vaccination in some cases. SARS-CoV-2 antibody specificity, breadth and 
maturation are affected by imprinting from exposure history, and distinct histological and 


antigenic contexts in infection compared to vaccination. 
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Introduction 

The urgent need for countermeasures against the coronavirus disease 2019 (COVID-19) pandemic 
has spurred rapid development of SARS-CoV-2 vaccines of diverse formulations. mRNA vaccines 
BNT162b2 (BioNTech-Pfizer) and mRNA-1273 (Moderna/NIAID) have demonstrated high 
efficacy and safety in clinical trials for COVID-19 prevention (Baden et al., 2021; Polack et al., 
2020; Walsh et al., 2020). Additional COVID-19 vaccines including adenoviral vectored vaccines 
ChAdOx1-S (AstraZeneca) (Voysey et al., 2021), Ad26.COV2.S (Johnson & Johnson) (Sadoff et 
al., 2021), and Gam-COVID-Vac (Sputnik V), and inactivated viral vaccines such as BBIBP-CorV 
(Sinopharm), also have reported efficacy. Correlates of vaccine-elicited protection from COVID- 
19 are the titers of neutralizing antibodies to SARS-CoV-2, and the concentration of antibodies 
binding to spike or receptor binding domain (RBD) (Earle et al., 2021; Gilbert et al., 2022; Khoury 
et al., 2021; Róltgen and Boyd, 2021). Most neutralizing antibodies target the RBD and prevent 
binding to angiotensin-converting enzyme 2 (ACE2) receptor (Greaney et al., 2021a; Yuan et al., 
2021). Current SARS-CoV-2 vaccines all contain or induce expression of antigens similar to those 
of the early Wuhan-Hu-1 viral isolate, but differ in elicited binding and neutralizing antibody 
responses, with higher responses from mRNA vaccines compared to adenovirus-vectored or 
inactivated virus vaccines (Dashdorj et al., 2021a, 2021b). It remains to be determined precisely 
how the immune system responds to mRNA and other vaccine platforms compared to SARS-CoV- 
2 infection. Data from RBD variant antigen yeast display and pseudotyped virus neutralization 
show that RBD epitope targeting by polyclonal serum antibodies is narrower in infected patients 
compared to mRNA-1273 vaccinees (Greaney et al., 2021b). 

Several SARS-CoV-2 variants of concern with mutations in the spike gene have emerged and 


spread globally, with differing abilities to evade neutralizing antibody responses elicited by 
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Wuhan-Hu-1 infection or vaccination. The most immune-evasive variants, including the recent 
Omicron variant, have alterations in epitopes containing amino acid E484 (Garcia-Beltran et al., 
2021; Greaney et al., 2021a; Hoffmann et al., 2021). The appearance of virus variants, waning 
antibody levels after infection or vaccination (Falsey et al., 2021; Levin et al., 2021), and 
breakthrough infections in previously immunized individuals (Keehner et al., 2021) indicate that 
periodic vaccine boosting of immunity to SARS-CoV-2 is warranted. Third doses of mRNA-1273 
(Chu et al., 2021) and BNT162b2 (Falsey et al., 2021) administered several months after the 
second dose prompt an increase in neutralizing antibodies greater than the peak following initial 
vaccination doses. mRNA-1273 vaccination followed by mRNA booster vaccines expressing Beta 
spike gives higher neutralizing titers to Wuhan-Hu-1-like SARS-CoV-2 compared to the Beta 
variant (Choi et al., 2021; Wu et al., 2021), suggesting that some degree of immune imprinting, or 
preferential responses to the viral variants initially encountered by the immune system, may affect 


the development of antibodies against new viral variants (Wheatley et al., 2021). 


Germinal center (GC) responses in human lymphoid tissues enable antibody affinity maturation 
and durable serological and memory B cell responses, although extrafollicular B cell responses are 
also reported (Elsner and Shlomchik, 2020; Lam et al., 2020; Woodruff et al., 2020). The degree 
to which SARS-CoV 2 infections or different vaccines stimulate GC responses and differ in factors 
such as the quantity, persistence and localization of antigen in lymph nodes (LNs) and other 
lymphoid tissues are important open questions. Approaches such as fine-needle aspiration (FNA) 
are being increasingly used to study LN-derived cells from healthy human subjects (Havenar- 
Daughton et al., 2020; Lederer et al., 2021; Turner et al., 2021). Disrupted LN GCs have been 


reported in autopsies of deceased COVID-19 patients (Haslbauer et al., 2021; Kaneko et al., 2020), 
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while elevated frequencies of GC B cells are seen after mRNA vaccination in healthy individuals 
(Turner et aL, 2021) and lower GC B cell frequencies after mRNA vaccination of 
immunocompromised individuals (Lederer et al., 2021). To date, no direct comparison of LN GC 
histology and cellular composition, combined with measurement of viral or vaccine antigen 
quantity, persistence and distribution in draining LN sites of COVID-19 patients and vaccinees 


has been reported. 


Неге, we compare antibody responses in BNT162b2 mRNA vaccine recipients following 1“, 214 
and 314 vaccine doses, to antibody responses of COVID-19 patients. We find differences in the 
magnitude, isotype profiles, SARS-CoV-2 spike domain specificity and breadth of binding 
antibody responses to a panel of nine viral variants in addition to Wuhan-Hu-1. Anti-RBD IgG 
binding to SARS-CoV-2 variants of concern and interest in recipients of four different vaccines 
(BNT162b2, ChAdOx1-S, Gam-COVID-Vac, and BBIBP-CorV) and іп COVID-19 patients, 
shows greater binding breadth for viral variant RBDs following all vaccines compared to Wuhan- 
Hu-1 infection. We quantify a strong imprinting effect of prior vaccination with Wuhan-Hu-1 
spike antigen on antibody specificities following breakthrough infection with viral variants. 
Histological analysis of draining LN shows marked impairment of GCs in severe COVID-19 
compared to mRNA vaccination, higher quantities and persistence of spike antigen accumulated 
in the GCs of mRNA vaccinees, and detectable vaccine RNA in GCs for up to two months post- 


2™ dose. 


Results 
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Magnitude and waning of anti-SARS-CoV-2 IgG following BNT162b2 vaccination, and 
response to 374 dose boost 

We measured anti-SARS-CoV-2 antibodies for nucleocapsid (N), full spike and КВР in Stanford 
BNT162b2 study participant plasma samples using multiplexed electrochemiluminescence (ECL) 
assays (Meso Scale Discovery, MSD), in WHO Binding Antibody Units (BAU). The first and 
second vaccine doses were at day 0 and day 21, with third dose boosting at approximately 9 
months. Plasma samples were collected in a time course up to seven months after the first dose, 
and up to one month after the third dose. Four of the 59 vaccine recipients had a history of SARS- 
CoV-2 reverse-transcription quantitative polymerase chain reaction (RT-qPCR)-confirmed 
infection (CoV-2+) prior to vaccination. IgG for spike protein and КВР in vaccinees reached their 
initial peak at day 28 after the first dose (Figures 1A and 1B). IgG binding to spike and RBD was 
highly correlated with SARS-CoV-2 neutralization titers (Arunachalam et al., 2021) (Figure ІС). 
By 9 months, spike-specific IgG had decreased approximately 20-fold from the maximum, but the 
3" dose boost raised IgG concentrations above the prior peak within one week. IgG specific for N 
protein, which is not encoded in the vaccine, was negative throughout the study in 54 of the 
previously uninfected vaccinees, but one participant seroconverted for anti-N IgG between day 90 
and 210 after the prime, indicating a breakthrough infection (Figure 1A). СоУ-2-- vaccinees had 
accelerated RBD and spike IgG responses after the first dose, and detectable anti-N IgG unaffected 


by vaccination (Figure S1A). 


BNT162b2 recipients had weak IgM and IgA responses to spike and RBD compared to their IgG 
responses. Robust IgG responses were seen in all age groups (Figure 1B; Figures 51А - SIC). 


Convalescent COVID-19 patients and BNT162b2 vaccinees had similar low saliva IgG 
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concentrations for spike апа КВР, several orders of magnitude lower than those detected in plasma 
(Figure S1D). As in plasma, saliva IgG peaked at one week after 34 dose boosting, at higher values 
than the peak after the 274 dose (Figure S1D). Reported side-effects after vaccination showed no 


relationship to plasma IgG responses (Figures 52А and S2B). 


BNT162b2 vaccination and Wuhan-Hu-1 SARS-CoV-2 infection stimulate distinct antibody 
isotypes and endemic coronavirus antibody responses 

Severe COVID-19 stimulates higher SARS-CoV-2-specific antibody titers than asymptomatic 
infection or mild illness (Long et al., 2020; Róltgen et al., 2020). We compared antibody isotype 
concentrations specific for spike and КВР” in COVID-19 patients (Stanford cohort 1 of this study) 
from the initial months of the pandemic (Róltgen et al., 2020), to the responses of the Stanford 
BNT162b2 vaccinees (Figures 2A and 2B). Patients were classified as outpatients; admitted 
patients not requiring care in the intensive care unit (ICU); ICU patients; and those who died from 
their illness. Stanford BNT162b2 vaccinee RBD and spike IgG concentrations were comparable 
to those of severely ill patients, and higher than those of mildly or moderately ill patients for anti- 
RBD antibodies at day 42 (Figures 2A and 2B). The BNT162b2 vaccine induced a highly IgG- 
polarized serological response with minimal IgM and IgA binding spike and RBD (Figures 2A 
and 2B). Principal component analysis (PCA) showed clustered and homogeneous SARS-CoV-2 
spike and spike domain-specific serological responses in BNT162b2 vaccinees compared to 
infected patients, as quantified by smaller distances for vaccinated participants from the group 
centroid (Figures 2C and 2D). 

SARS-CoV-2 vaccinees and COVID-19 patients showed boosting of SARS-CoV-1 spike 


antibodies, but infected patients showed greater boosting of spike IgG and IgA for endemic human 
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betacoronaviruses OC43 and HKUI (Figures S3A апа S3B). Тһе 34 BNT162b2 vaccine dose 
further increased vaccinee titers to SARS-CoV-1, ОСАЗ, and HKUI (Figure S3C). Antibodies to 
the spike antigens of the endemic human alphacoronaviruses NL63 and 229E were not boosted 


(Figure S3). 


Greater breadth of IgG binding to viral variants following BNT162b2 vaccination compared to 
infection with Wuhan-Hu-1 SARS-CoV-2 

Immune-evasive SARS-CoV-2 variants have spread globally (Harvey et al., 2021; Plante et al., 
2021; Róltgen and Boyd, 2021). We compared plasma IgG responses to the RBDs of nine different 
SARS-CoV-2 variants of concern and interest in BNT162b2 vaccinees and COVID-19 patients, 
using multiplexed MSD ECL assays. For RBD antigens from Epsilon, Kappa, B.1.526.2, 
B.1.214.2, Alpha, Eta/Iota, Gamma, P.3, and Beta variants, both vaccinee and infected patient IgG 
showed the greatest decrease in binding to Beta, Gamma, and P.3 variants relative to Wuhan-Hu- 
1 (Figure 3A). To quantify the differences in variant RBD binding by vaccinee and patient plasma 
IgG we calculated the ratios of anti-RBD IgG concentrations for Wuhan-Hu-1 compared to viral 
variants, with higher ratios indicating greater binding of Wuhan-Hu-1 RBD compared to variant 
RBD (Figure 3B). COVID-19 patients showed a greater IgG binding bias for Wuhan-Hu-1 RBD 
compared to variant RBDs in the initial weeks post-onset of symptoms; in contrast, BNT162b2 
vaccinee IgG had relatively greater breadth of binding to variant RBDs and less preference for 
Wuhan-Hu-1 RBD. Over time, infected patient plasma samples showed improvement in variant 
RBD binding relative to Wuhan-Hu-1 RBD, suggesting evolution of the antibody response through 
at least 7 weeks post-onset of symptoms (Figure 3). BNT162b2 vaccinee IgG Wuhan-Hu-1 to 


variant RBD binding ratios did not change from day 21 onward. The greater breadth of variant 
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RBD binding (including the Delta variant) by vaccinee IgG compared to COVID-19 patient IgG 
was seen in a second, independent validation cohort (Stanford cohort 2) of predominantly mildly 
ill COVID-19 patients. Greater Wuhan-Hu-1 to variant RBD IgG binding ratios were found in 
week 2 to 3, month 1, month 3, and month 7 in infected patients compared to vaccinees (Figures 
S4A and S4B), with improvement in variant recognition over time in the infected patients. Notably, 
the increased breadth of vaccinee IgG compared to COVID-19 patient IgG binding to viral variant 
antigens was greatest for RBD, the main target of neutralizing antibodies, and was decreased or 
not detected when whole spike antigens were tested (Figure S4C). Functional blocking of ACE2 
binding to RBD was concordant with the RBD-specific IgG concentrations measured in these 


populations (Figures S4A and S4D). 


Improved IgG binding to viral variants is consistent across four COVID-19 vaccines 
(BTNI62b2, ChAdOx1-S, Gam-COVID-Vac and BBIBP-CorV) compared to infection 

Several COVID-19 vaccines, including mRNA, viral vector-based, and inactivated virus vaccines, 
have been approved for use internationally. Varying efficacy and antibody responses from the 
vaccines have been reported (Baden et al., 2021; Dashdorj et al., 2021a, 2021b; Polack et al., 2020; 
Sadoff et al., 2021; Voysey et al., 2021). We compared IgG responses in Stanford COVID-19 
cohort 2 patients and BNT162b2 vaccinees to those of participants in a Mongolian observational 
study deploying four different COVID-19 vaccines: the mRNA vaccine BNT162b2; adenoviral 
vectored vaccines ChAdOx1-S (AstraZeneca) and Gam-COVID-Vac (Sputnik V); and an alum- 
adjuvanted, inactivated viral vaccine BBIBP-CorV  (Sinopharm). RBD-specific IgG 
concentrations for Wuhan-Hu-1 and all viral variants measured (Epsilon, Kappa, B.1.526.2, Delta, 


Alpha, Eta/Iota, Gamma, P.3, and Beta) differed greatly between vaccine groups, with BNT162b2 
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eliciting the highest antibody levels, followed by AstraZeneca, Sputnik V, and Sinopharm 
vaccination (Figure 4A). IgG concentration differences between vaccines were significant for most 
viral variant RBDs. Stanford BNT162b2 vaccinees compared to Mongolian BNT162b2 vaccinees 
had higher IgG concentrations at early time points, likely due to differences in timing of sample 
collection (Stanford day 28 and day 90; Mongolian participants variable time points before 1- 
month and 3-months) (Figure 4A). Despite the different vaccine compositions and magnitudes of 
antibody responses, all four vaccines elicited IgG with relatively greater breadth of viral variant 


RBD binding compared to that of infected patients (Figure 4B). 


Variant-specific serological responses following Alpha and Delta SARS-CoV-2 infection, and 
immune imprinting after vaccination 

Immune imprinting, a phenomenon in which primary exposure to an antigen forms epitope- 
specific B cell memory and affects future B cell and antibody responses against variant epitopes, 
has been studied in influenza infection and vaccination. COVID-19 patients and BNT162b2 
vaccinees who were only exposed to Wuhan-Hu-1 antigens in this study exhibit a consistent 
hierarchy in IgG binding concentrations to the different SARS-CoV-2 variant RBDs relative to the 
Wuhan-Hu-1 RBD, decreasing from Epsilon, Kappa, B.1.526.2, Delta, Alpha, Eta/Iota, Gamma, 
P.3, to Beta (Figure S5A). To test for imprinting of the serological response to variant RBDs, we 
first analyzed the ratios of Wuhan-Hu-1 to variant RBD IgG concentrations in COVID-19 patients 
who were infected with Alpha or Delta variants, confirmed by allele-specific RT-qPCR testing or 
viral sequencing. IgG from Alpha or Delta variant-infected patients with no history of COVID-19 
vaccination or prior SARS-CoV -2 infection preferentially bound Alpha and Delta variant RBDs, 


respectively, compared to Wuhan-Hu-1 RBD (Figure 5A, upper panels). Delta infection also 
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elicited higher IgG concentrations to other variant RBDs containing L452R such as Epsilon and 
Kappa, compared to Wuhan-Hu-1 (Figure 5A, upper right). PCA of variant RBD-specific IgG 
responses in vaccinees and variant-infected patients (Figure S5B) highlights the distinct 
serological responses elicited by infection with the variant viruses. To test whether prior exposure 
to one SARS-CoV-2 RBD variant causes imprinting of humoral immunity, we analyzed plasma 
from individuals vaccinated with Wuhan-Hu-1-like antigens and subsequently infected with Alpha 
or Delta variants (Figure 5A, lower panels). Despite breakthrough infection with Alpha or Delta 
viral variants, the vaccinated individuals showed patterns of IgG binding to viral variant RBDs 
similar to those of individuals exposed to only Wuhan-Hu-1. We quantified the degree of 
imprinting of IgG specificity by log-transforming the ratios of IgG binding to pairs of antigens (for 
example, Wuhan-Hu-1 RBD compared to Delta RBD) for individual samples, then rescaling to 
range from -100% to +100% corresponding to the maximal preference for each antigen observed 
in other plasma specimens, including those from individuals exposed only to a single antigen 


variant (Figure 5B). 


LN GC impairment in severe SARS-CoV-2 infection but robust development following SARS- 
CoV-2 mRNA vaccination 

The differences in viral variant RBD IgG binding between SARS-CoV-2 infected patients and 
recipients of the four COVID-19 vaccines suggest that the organization of the humoral immune 
responses in secondary lymphoid tissues may differ between infection and vaccination, potentially 
due to direct effects of the viral infection, differences in innate immune stimuli between 
vaccination and infection, or the quantity or localization of viral antigens, among other 


possibilities. Previous studies have revealed a loss of GCs and a reduction in BCL6* GC B cells 
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in severe acute SARS-CoV-2 infection, raising the possibility that humoral responses may be 
altered or subverted by the virus (Kaneko et al., 2020). It is unclear whether draining LN immune 
responses to SARS-CoV-2 infection in the lungs differ from those elicited in axillary LNs 
following deltoid intramuscular mRNA vaccination. To compare GC architecture in response to 
SARS-CoV-2 infection and vaccination, we obtained peribronchial LN tissues from six COVID- 
19 patients and three control autopsy cases as well as axillary LN core needle biopsies of seven 
individuals vaccinated with mRNA-1273 or BNT162b2. Importantly, core needle biopsy sample 
tissue volumes were suitable for assessment of LN histoarchitecture. Vaccinee axillary LN 
biopsies were from the ipsilateral (same-side) arm vaccinated. Controls were thoracic LNs from 
individuals who succumbed to pre-pandemic non-COVID-19 pneumonias, and contain GCs likely 
due to ongoing adaptive immune responses elicited by other antigens. LN histology for COVID- 
19 patients and vaccinees was evaluated with 4-color co-detection by indexing (CODEX) 
immunofluorescence analysis for CD20, CD3, BCL6 and CD21, which are markers of B cells, T 
cells, GC B cells (or T follicular helper (Tfh) cells) and follicular dendritic cells, respectively 
(Figure 6A, Figure S6A), as well as by single-color immunohistochemical stains for these markers 
and (һе Tfh cell marker PD-1 (Figure 6B). GCs were poorly formed in the severely ill COVID-19 
patient peribronchial LNs compared to the axillary LNs of vaccinees, with disrupted CD21+ 
follicular dendritic cell networks and decreased BCL6* cells (including GC B cells and Tfh cells) 
and PD-1* cells (consistent with Tfh cells) within GCs (Figures 6A - 6E). Disruption of CD21* 
follicular dendritic cell networks was seen in both primary and secondary follicles in LNs from 
COVID-19 patients (Figure S6B). mRNA vaccination was associated with follicular hyperplasia 
with fully developed GC architecture, including robust induction of GC B cells, Tfh cells and 


extensive follicular dendritic cell networks (Figures 6A and 6B). 
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Prolonged detection of vaccine mRNA іп LN GCs, and spike antigen in LN GCs and blood 
following SARS-CoV-2 mRNA vaccination 

The biodistribution, quantity and persistence of vaccine mRNA and spike antigen after 
vaccination, and viral antigens after SARS-CoV 2 infection, are incompletely understood but are 
likely to be major determinants of immune responses. We performed in situ hybridization with 
control and SARS-CoV-2 vaccine mRNA-specific RNAScope probes in the core needle biopsies 
of the ipsilateral axillary LNs that were collected 7-60 days after 2"! dose of mRNA-1273 or 
BNT162b2 vaccination, and detected vaccine mRNA collected in the GCs of LNs on day 7, 16, 
and 37 post vaccination, with lower but still appreciable specific signal at day 60 (Figures 7A - 
7E). Only rare foci of vaccine mRNA were seen outside of GCs. Axillary LN core needle biopsies 
of non-vaccinees (n - 3) and COVID-19 patient specimens were negative for vaccine probe 
hybridization. Immunohistochemical staining for spike antigen in mRNA vaccinated patient LNs 
varied between individuals, but showed abundant spike protein іп GCs 16 days post-2"4 dose, with 
spike antigen still present as late as 60 days post-2"! dose. Spike antigen localized in a reticular 
pattern around the GC cells, similar to staining for follicular dendritic cell processes (Figure 7B). 
COVID-19 patient LNs showed lower quantities of spike antigen, but a rare GC had positive 
staining (Figure 7F). Immunohistochemical staining for N antigen in peribronchial LN secondary 
and primary follicles of COVID-19 patients (Figures 7F - 7I) was positive in 5 of the 7 patients, 


with a mean percentage of nucleocapsid-positive follicles of more than 25%. 


Spike protein was detected in the plasma of 96% of the vaccinees at days 1-2 (median spike 


concentration of 47 pg/mL) and in 63% at day 7 (median spike concentration of 1.7 pg/mL) after 
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the prime vaccine dose. In contrast, spike antigen detection after the vaccine boost on day 21 was 
reduced, with half of the study participants being positive on day 1-2 (median spike concentration 
of 1.2 pg/mL) and only one individual on day 7 post-boost (Figure 7J). We suspected that high 
concentrations of spike-specific antibodies developed by vaccinees within the first 2 to 3 weeks 
after the prime vaccine dose could impede detection of spike antigen by competing for spike 
binding sites with the anti-spike reagent antibodies in the antigen assay. To test this hypothesis, 
we added different concentrations of recombinant spike protein to spike-negative vaccinee plasma 
samples collected on day 0 (n = 3) and day 28 (n = 3) after the prime vaccination. While the 
recombinant spike protein could readily be detected in day O plasma samples, only high 
concentrations of the antigen led to a positive signal when mixed with the day 28 samples (Figure 
7K). We then mixed a spike-positive plasma sample collected опе day after vaccination with spike- 
negative plasma samples collected on days 0, 21, 22-23, and 28 (n - 4 each). Spike antigen 
detection levels were high in the mix of day 1 and day 0 samples, decreased in the mix of day 1 
and day 21, and day 1 and day 22-23 samples, and below the cutoff for positive in the mix of day 
] and day 28 samples (Figure 7L). Together our results are consistent with spike-specific 


antibodies blocking detection of the antigen in antigen capture-based assays. 


Discussion 

One of the positive developments amid the global calamity of the SARS-CoV-2 pandemic has 
been the rapid design, production and deployment of a variety of vaccines, including remarkably 
effective mRNA vaccines encoding the viral spike (Baden et al., 2021; Polack et al., 2020). We 
find that BNT162b2 vaccination produces IgG responses to spike and RBD at concentrations as 


high as those of severely ill COVID-19 patients and follows a similar time course. Unlike infection, 
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which stimulates robust but short-lived IgM and IgA responses, vaccination shows a pronounced 
bias for IgG production even at early time points. These responses were similar across the adult 
age range in our study. The relative absence of IgM and IgA responses suggests a potent effect of 
the vaccine formulation in driving early and extensive IgG class-switching, potentially as a result 
of the reported T helper type 1-polarized CD4* T cell responses stimulated by vaccine components 
(Lederer et al., 2020; Lindgren et al., 2017; Pardi et al., 2018). Our data demonstrate that vaccinee 
plasma and saliva spike and RBD-specific IgG concentrations decrease from their peak values by 
approximately 20-fold by 9 months after primary vaccination, but quickly exceed prior peak 


concentrations in 7 to 8 days after boosting with a 3 vaccine dose. 


Correlates of immunological protection from SARS-CoV-2 infection following vaccination or 
prior infection are still under investigation. Analysis of Moderna mRNA-1273 and AstraZeneca 
ChAdOx1-S responses highlights the overall similarity of correlate of protection results for spike- 
binding antibody and neutralizing antibody assays (Feng et al., 2021; Gilbert et al., 2022). We 
compared spike or RBD-binding antibody responses to Wuhan-Hu-1 SARS-CoV-2 neutralization 
data in BNT162b2 vaccinees and confirmed the high correlation of these assay results, supporting 
the interpretation that sensitive, precise and validated commercial multiplexed antigen-binding 
assays with a wide dynamic range, such as the MSD ECL assays in this study, will be valuable in 
providing standardized correlates of protection data for vaccines as the pandemic continues. 
Particularly in the context of viral variants, it will be important to determine whether predictions 
of vulnerability to infection or severe disease can be improved by adding data from other 


immunological assays, including T cell measurements. 
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Differences in B cell responses to SARS-CoV-2 infection and vaccination may be reflected in the 
binding breadth of antibodies to different SARS-CoV-2 variants. We find that plasma of 
individuals who received prime/boost BNT162b2 vaccination, as well as individuals who received 
adenoviral vectored (ChAdOx1-S or Gam-COVID-Vac) or inactivated virus (BBIBP-CorV) 
COVID-19 vaccines show consistent patterns of binding to variant RBDs with modest decreases 
compared to Wuhan-Hu-1 RBD binding. In contrast, COVID-19 patients produce antibody 
responses with significantly greater Wuhan-Hu-1 RBD binding preference and lower breadth of 
variant RBD binding. These differences between vaccinee and COVID-19 patient IgG variant 
antigen binding were greatest for the RBD, the target of most neutralizing antibodies, and were 
diminished when full spike antigen with its greater number of non-neutralizing epitopes was tested. 
These results, covering many clinically relevant viral variant antigens and several vaccine 
modalities, are consistent with findings for RBD binding IgG in mRNA-1273 vaccinees compared 
to infected patients (Greaney et al., 2021b). Notably, COVID-19 patients with Alpha or Delta 
variant infections display characteristic serological profiles specific to the RBD of the infecting 
variant, indicating that SARS-CoV-2 variant serotyping may be useful for epidemiological studies 
of populations to determine exposure to circulating SARS-CoV-2 variants. Both vaccinees and 
COVID-19 patients exposed to Wuhan-Hu-1 antigens show the greatest decreases in antibody 
binding to RBD variants harboring E484 alterations, including Beta and Gamma. Although 
susceptibility to infection by viral variants is common to both vaccinated and convalescent 
populations, particularly as antibody titers decrease over time (Israel et al., 2021; Levin et al., 
2021), our findings lead to the prediction that antibodies derived from infection may provide 
somewhat decreased protection against virus variants compared to comparable concentrations of 


antibodies stimulated by vaccination. 
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As additional variants of SARS-CoV-2 appear over time, individuals will acquire distinct 
immunological histories depending on which vaccines they received and which viral variants 
infected them. The idea that “imprinting” by a prior antigen exposure can shape, either positively 
or negatively, the response to a subsequent variant is well established in studies of influenza 
viruses, and has been implicated in birth-year differences in susceptibility to particular avian 
influenza viruses (Gostic et al., 2016). We find that prior vaccination with Wuhan-Hu-1-like 
antigens followed by infection with Alpha or Delta variants gives rise to plasma antibody 
responses with apparent Wuhan-Hu-1-specific imprinting manifesting as relatively decreased 
responses to the variant virus epitopes, compared to unvaccinated patients infected with those 
variant viruses. While current booster vaccinations are still based on the Wuhan-Hu-1-like 
antigens, vaccine manufacturers are in the process of evaluating updated vaccines encoding 
sequences from one or more circulating variants. Initial results from 214 dose boosting with Beta 
spike-encoding mRNA vaccines after prior 2-dose mRNA-1273 vaccination are consistent with 
our findings of significant imprinting of serological responses by the first antigen encountered 
(Choi et al., 2021; Chu et al., 2021), indicating that vaccine-derived imprinting affects subsequent 
antibody responses stimulated by vaccination as well as infection. The extent to which vaccine 
boosting or infection with different variants will effectively elicit antibody responses to new 
epitopes, or rather increase responses to the epitopes of antigens encountered previously, as in the 
"original antigenic sin" phenomenon described for influenza virus infection and vaccination 
(Arevalo et al., 2020; Zhang et al., 2019), will be an important topic of ongoing study. The degree 
of imprinting may depend on the particular variants and the order in which they are introduced to 


the individual's immune system, and the number of exposures, such as the number of vaccine 
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doses received. Additional data for evaluating the magnitude of these effects and their 
consequences for protection from infection are likely to become available in coming months, as 
individuals with different histories of SARS-CoV-2 vaccination or viral variant infection become 
infected with the more highly mutated Omicron variant (https://covdb.stanford.edu/page/mutation- 
viewer/#omicron). As a practical consideration, the very high spike-specific IgG concentrations 
generated by mRNA vaccination and periodic additional booster doses may be able to compensate 
for relatively decreased binding to new viral variant antigens, potentially decreasing the public 


health impact of antibody response imprinting if vaccine boosting is widely adopted. 


We hypothesized that differences in the serological responses observed in SARS-CoV-2 infection 
compared to vaccination, particularly those related to variant antigen binding breadth, could be 
related to the anatomical sites where the viral antigens are encountered, the quantity of viral 
antigens in those anatomical sites, differences in the cell populations stimulated in secondary 
lymphoid tissues, and potential damage to immunological tissues during infection. With CODEX 
multiplexed immunofluorescence microscopy and immunohistochemical microscopy, we 
identified follicular hyperplasia with robust axillary LN GCs after mRNA (BNT162b2 or mRNA- 
1273) vaccination, containing CD21- follicular dendritic cell networks, BCL6* B cells and PD-1* 
cells at significantly higher frequencies compared to those in peribronchial LNs of deceased 
COVID-19 patients. These findings demonstrate greater stimulation of GC B cells and Tfh cells 
in vaccination, and normal functional organization of GC follicular dendritic cells. Loss or 
impairment of GCs in patients with severe COVID-19 suggests that SARS-CoV-2 viral infection 
subverts the humoral immune response, by directly damaging immune cells or as a secondary 


effect of inflammatory responses to infection (Feng et al., 2020; Kaneko et al., 2020). The observed 
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extended presence of vaccine mRNA and spike protein in vaccinee LN GCs for up to 2 months 
after vaccination was in contrast to rare foci of viral spike protein in COVID-19 patient LNs. We 
hypothesize that the abundant spike antigen in the GCs of mRNA vaccine recipient LNs may 
contribute to the increased breadth of viral variant RBD binding by IgG seen after vaccination, 
potentially due to high antigen concentrations stimulating B cells with lower affinity for Wuhan- 
Hu-1 spike epitopes and better binding to variant epitopes. Persistent vaccine RNA and spike 
antigen at elevated concentrations in vaccinee LNs could result in less strict selection for higher- 
affinity B cells in the immune response compared to situations where antigen is more limiting 
(Cirelli et al., 2019). However, our observation that all vaccine modalities (mRNA, adenoviral and 
inactivated virus) stimulated greater viral variant breadth of IgG binding than infection could 
indicate that some other aspect of SARS-CoV-2 infection underlies these differences, such as 


alteration of GC function. 


Pre-pandemic analysis of a model RNA vaccine for yellow fever virus in a rhesus macaque at 16 
hours post-vaccination showed that vaccine RNA іп LN cell suspensions was detected 
predominantly in professional antigen-presenting cells including monocytes, classical dendritic 
cells and B cells at this early time point (Lindsay et al., 2019). Data from follicular dendritic cells 
were not reported. Our histological data from SARS-CoV-2 mRNA-vaccinated humans at 
considerably later time points (7 to 60 days post-2™ dose) show vaccine RNA almost entirely in 
GCs, distributed primarily between the nuclei of GC cells, similar to the pattern seen by 
immunostaining for follicular dendritic cell processes or B cell cytoplasm. Additional co- 
localization studies with higher resolution may be required to determine more exactly which 


specific cell types harbor mRNA vaccine and spike antigen in humans following COVID-19 
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mRNA vaccination and infection, and may provide further mechanistic insights into the basis for 


the differences in serological responses after vaccination compared to infection. 


At least some portion of spike antigen generated after administration of BNT162b2 becomes 
distributed into the blood. We detected spike antigen in 9646 of vaccinees in plasma collected one 
to two days after the prime injection, with antigen levels reaching as high as 174 pg/mL. The range 
of spike antigen concentrations in the blood of vaccinees at this early time point largely overlaps 
with the range of spike antigen concentrations reported in plasma in a study of acute infection 
(Ogata et al., 2020), although a small number of infected individuals had higher concentrations in 
the ng/mL range. At later time points after vaccination, the concentrations of spike antigen in blood 
quickly decrease, although spike is still detectable in plasma in 6396 of vaccinees one week after 
the first dose. A practical finding in our study is that the detection of spike antigen in plasma 
samples is impeded after 2"! dose BNT162b2 vaccination, likely due to the formation of circulating 
immune complexes of anti-spike antibodies and spike protein, masking the antigen epitopes of the 
capture and detection antibodies that form the basis of antigen detection assays, similar to assay 
interference that has been reported for other diseases (Bollinger et al., 1992; Lima et al., 2014; 


Miles et al., 1993). 


Limitations of the study 

Data from SARS-CoV-2 infected clinical cohorts and vaccinated individuals in this study are 
observational. Longitudinal data for COVID-19 vaccine responses are derived predominantly from 
BNT162b2 mRNA vaccine recipients at Stanford, with data for the other four COVID-19 vaccines 


at a single post-vaccination time point per individual. To make precise, internally controlled 
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comparisons of polyclonal antibody responses to different viral variant antigens, we used 
multiplexed ECL assays of antibody binding to RBD, rather than virus neutralization assays, 
therefore our data do not reflect potentially functional antibodies binding to the spike N-terminal 
domain, or antibodies that may have other activities in vivo. In the analysis of imprinting of 
serological responses, plasma specimens were not available from the period after vaccination but 
before variant virus infection, precluding direct comparison of antibody specificities pre- and post- 
infection. Additional epidemiological studies will be needed to evaluate the clinical impact of 
antibody response imprinting on susceptibility to infection by new viral variants, and the severity 
of disease in infected patients. LN histological comparisons between COVID-19 patients and 
vaccinees have the limitations that the infected patient specimens were limited to those with severe 
disease; the number of individuals analyzed was relatively low (six COVID-19 patients and seven 
vaccinees); and the LN sampling was not done prospectively at pre-determined time points after 
vaccination or infection. The serological analysis in this study is of polyclonal antibody responses; 
evaluation of the clonal B cell and plasma cell populations producing these antibody mixtures in 
comparable numbers of subjects in infection and vaccination will likely be required for further 


mechanistic understanding. 


Taken together, these results underscore important differences between SARS-CoV-2 antibody 
responses produced by vaccination versus infection. Key questions for the months and years ahead 
include the duration of effective vaccine-stimulated serological responses after 3 dose boosting 
or other repeated exposures, particularly for the recent Omicron variant and other variants that will 
emerge in future, and the safety and efficacy of variant-targeting vaccine boosters in previously 


vaccinated or infected individuals. Further mechanistic investigations into the differences in 
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antibody breadth elicited by vaccination and infection are needed to define the roles of T cell help, 
antibody affinity maturation, GC function, and innate immune responses to vaccine components, 
as well as the cellular and subcellular distribution of vaccine RNA and expressed antigen in 
Iymphoid tissues. Lessons from the antibody responses to the initial SARS-CoV-2 variants are 
likely to be important both for preparing for future additional variants of this virus, as well as 


improving vaccination strategies for other pathogens such as influenza virus. 
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Figure titles and legends 
Figure 1. Magnitude and duration of anti-SARS-CoV-2 IgG following BNT162b2 
vaccination and 3"! dose boost. 

(A) Anti-SARS-CoV-2 М, КВР, and spike (S) antibody responses are shown for plasma samples 
from individuals who received BNT162b2 prime (DO, n = 59 individuals), second dose (D21, n = 
58 individuals), and third dose (around month 9, n - 36 individuals) vaccination. Box-whisker 
plots of the anti-SARS-CoV-2 IgG concentrations in WHO BAU/mL show the interquartile range 
as the box and the whisker ends as the most extreme values within 1.5 times the interquartile range 
below the 25% quantile and above the 75% quantile. Red dashed lines indicate the cutoff values 
for positivity of each assay (MSD, package insert). 

(B) Heatmap showing the development of antibody responses in longitudinal samples collected at 
DO, D7, D21, D28, D42, and D90/120 time points post-prime vaccination (x-axis). WHO BAU/mL 
Ig concentrations are displayed for study participants sorted by age (y-axis, color-coded). Rows 
are labeled on the right with “СоУ-2--” for participants with a previous SARS-CoV-2 RT-qPCR 
positive test result. 

(C) Correlations between anti-RBD and anti-spike IgG binding antibody concentrations in WHO 
BAU/mL and SARS-CoV-2 virus neutralization assays. Spearman rank correlation (coefficient = 
Rho, displayed in the plot for each assay comparison) was used to assess the strength of correlation 
between binding antibody concentrations and virus neutralization results. Red dashed lines 


indicate the cutoff values for positivity of each assay (MSD, package insert). 


Figure 2. BNT162b2 vaccination and SARS-CoV-2 infection elicit distinct Ig isotype profiles. 


(A, B) Anti-SARS-CoV-2 N, RBD, and spike (S) IgM, IgG, and IgA antibody responses are shown 
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for individuals who received BNT162b2 prime (DO) and 274 (D21) vaccination doses and for 
COVID-19 patients. 

(A) The heatmap shows the development of antibody responses in longitudinal samples from 
vaccinees/patients collected at DO, D7 / week 1, D21 / weeks 2&3, D28 / week 4, D42 / weeks 
5&6, and D90/120 / >угеек 7 after vaccination / COVID-19 symptom onset (x-axis). The color 
scale encodes the median values of log10 WHO BAU/mL Ig concentrations. 

(B) Box-whisker plots show the development of antibody responses in longitudinal samples from 
vaccinees / patients collected at DO, D7 / week 1, D21 / weeks 2&3, D28 / week 4, D42 / weeks 
5&6, and D90/120 / >week 7 after vaccination / COVID-19 symptom onset (x-axis). Box-whisker 
plots show the interquartile range as the box and the whisker ends as the most extreme values 
within 1.5 times the interquartile range below the 25% quantile and above the 75% quantile. 
Statistical test: pairwise Wilcoxon rank sum test with Bonferroni correction. *p « 0.05, **p « 0.01, 
***p < 0.001. 

Individuals were classified as outpatients (Outpt) and hospital admitted patients (Admit); intensive 
care unit (ICU) patients and those who died from their illness (Death); and vaccinees who had 
(CoV-2+) or had not had a positive SARS-CoV-2 test in the past. 

(C) PCA of anti-SARS-CoV-2 КВР, N-terminal domain, and S (but not М) IgM, IgG, and IgA 
concentrations across BNT162b2 vaccinees апа Wuhan-Hu-1-infected Stanford COVID-19 
patient cohort 1 at different time points after vaccination / COVID-19 symptom onset visualized 
on a consistent PCA reference created using D21 / weeks 2&3 as a reference time point. 

(D) Distribution of Euclidean distances between BNT162b2 vaccinee samples and their centroid, 
compared to Wuhan-Hu-l-infected Stanford COVID-19 patient cohort 1 samples and their 


centroid, at different time points after vaccination / COVID-19 symptom onset. 
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Figure 3. Greater breadth of IgG binding to SARS-CoV-2 RBD variants following 
BNT162b2 vaccination compared to infection with Wuhan-Hu-1 SARS-CoV-2. Anti-SARS- 
CoV-2 Wuhan-Hu-1 and viral variant RBD IgG responses are shown for Stanford individuals who 
received BNT162b2 vaccination and for Wuhan-Hu-1-infected COVID-19 Stanford patient cohort 
1 at different time points after vaccination / COVID-19 symptom onset. Box-whisker plots show 
the interquartile range as the box and the whisker ends as the most extreme values within 1.5 times 
the interquartile range below the 2596 quantile and above the 75% quantile. Significance between 
patient and vaccinee groups were tested with two-sided Wilcoxon rank sum test. *p < 0.05, **p « 
0.01, ***p « 0.001. 

(A) Anti-RBD IgG concentrations. 


(B) Ratios of anti-Wuhan-Hu-1 to variant КВР” IgG concentration. 


Figure 4. Greater breadth of IgG binding to SARS-CoV-2 variant RBDs following 
vaccination with four different vaccines compared to infection with Wuhan-Hu-1 SARS- 
CoV-2. Anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant RBD IgG responses are shown for 
individuals who received BNT162b2 (BioNTech-Pfizer), СҺАаОх1-5 (AstraZeneca), Gam- 
COVID-Vac (Sputnik V), and BBIBP-CorV (Sinopharm) vaccination and for Wuhan-Hu-1- 
infected COVID-19 Stanford patient cohort 2 within one month and around three months after 
vaccination / COVID-19 symptom onset. Box-whisker plots show the interquartile range as the 
box and the whisker ends as the most extreme values within 1.5 times the interquartile range below 
the 25% quantile and above the 75% quantile. Significance between groups were tested with 
pairwise Wilcoxon rank sum test with Bonferroni correction. *p « 0.05, **p « 0.01, ***p « 0.001. 


(A) Anti-RBD IgG concentrations. 
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(B) Ratios of anti-Wuhan-Hu-1 to variant КВР” IgG concentration. 


Figure 5. Variant-specific serological signature following Alpha and Delta SARS-CoV-2 
infection. 

(A) Anti-Wuhan-Hu-1 to variant RBD IgG concentration ratios are shown for individuals with 
primary SARS-CoV-2 Alpha or Delta variant infection (upper panels) or secondary variant 
infection after vaccination (lower panels). Box-whisker plots show the interquartile range as the 
box and the whisker ends as the most extreme values within 1.5 times the interquartile range below 
the 25% quantile and above the 75% quantile. 

(B) Anti-SARS-CoV-2 variant IgG binding preference levels of BNT162b2 vaccinees оп day 28 
post vaccination and of previously vaccinated or nonvaccinated individuals infected with the 


SARS-CoV-2 Delta variant. 


Figure 6. Disrupted LN GCs in COVID-19 patients versus mRNA vaccinees. 

(A) Representative LN GC histology of COVID-19 patients and vaccinees evaluated with 4-color 
CODEX immunofluorescence analysis for CD20 (B cells), CD3 (T cells), BCL6 (GC B cells 
(major subset) and follicular helper T cells (minor subset)), and CD21 (follicular dendritic cells). 

(B) Representative immunohistochemistry of GCs with CD21 (left), BCL6 (middle) and PD-1 
(right) in peribronchial LNs of an autopsy patient who died of COVID-19, a control autopsy patient 
who died from a non-COVID-19 pneumonia (pre-pandemic), and in an axillary LN of a patient 
vaccinated with a SARS-CoV-2 mRNA vaccine. 

(C-E) Relative proportion (upper) and absolute number (lower) of GCs in LNs (C), of BCL6+ cells 


within GCs (D), and of PD-1+ cells within GCs (E) from COVID-19 autopsy patients (n = 6), 


28 


631 


632 


633 


634 


635 


636 


637 


638 


639 


640 


641 


642 


643 


644 


645 


646 


647 


648 


649 


650 


651 


652 


control autopsy patients (n = 3), апа mRNA vaccinated patients (n = 7). Quantification performed 
in QuPath digital pathology analysis software. Wilcoxon rank sum test was used to calculate p 


values. Error bars represent mean + SEM. *p « 0.03; **p < 0.003. 


Figure 7. Localization of SARS-CoV-2 proteins and vaccine mRNA in LNs. 

(A) Representative LN GC after mRNA vaccination showing Hematoxylin & Eosin staining 
(upper left), 4-color CODEX staining (lower left), in situ hybridization of a SARS-CoV-2 mRNA 
vaccine-specific probe (upper right (lower magnification) and middle right (greater 
magnification)), and immunohistochemical (IHC) staining for spike antigen (lower right). Vaccine 
mRNA probe hybridization was visualized by colorimetric development with Fast Red 
chromogen, and positive IHC staining for spike antigen was visualized as granular brown color 
from 3,3'-Diaminobenzidine (DAB) reagent. 

(B) Representative in situ hybridization of an RNAScope control probe (left panels) and SARS- 
CoV-2 mRNA vaccine-specific probe (middle panels) within ipsilateral axillary core needle LN 
biopsies of female patients 7 to 60 days after second mRNA-1273 or BNT162b2 dose. Probe 
hybridization is indicated by red chromogen spots. IHC signal for spike antigen (right panels), is 
detected as granular brown staining. 

(C) Quantification of SARS-CoV-2 mRNA vaccine-specific probe-staining GCs in vaccinated LN 
biopsies. 

(D) Quantification of positive SARS-CoV-2 mRNA vaccine-specific probe spots per GC in 
vaccinee LNs. Error bars represent mean + SEM. 


(E) Spike protein-positive GC quantification from IHC staining of vaccinee LNs. 
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(F) IHC staining for spike (lower right panel) and nucleocapsid (upper panels and lower left panel) 
antigens in representative sections of COVID-19 patient peribronchial LNs. Nucleocapsid 
detection in primary (upper right panel) and secondary (upper left panel) LN follicles. 

(G) Due to the low frequency of detection of spike antigen in COVID-19 patient LNs, 
quantification is presented as the number of patients with positive staining in their LN specimens. 
(Н) Quantification of the number of COVID-19 patients with LN follicles positive for 
nucleocapsid IHC staining. 

(1) Number and percentage of nucleocapsid-positive follicles by IHC іп COVID-19 patient LNs. 
Error bars represent mean + SEM. 

(J) Spike concentration measured in plasma samples collected before and at several time points 
after BNT162b2 vaccination, with the red dotted line indicating the cutoff for positive. 

(K) Spike concentrations were measured in plasma samples collected from BNT162b2 vaccinees 
on DO (spike negative) or D28 (spike positive) spiked with different concentrations of recombinant 
spike protein. Black dotted line - cutoff for positive. 

(L) Spike concentration measured in plasma samples collected from BNT162b2 vaccinees on DO, 
D21, D22/23, and D28 mixed with the same plasma sample collected from one BNT162b2 


vaccinee on D1 (spike positive). Black dotted line = cutoff for positive. 
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STAR Methods 
RESOURCE AVAILABILITY 
Lead Contact 


Further information and requests for resources and reagents should be directed to the Lead Contact, 


Dr. Scott D. Boyd (publications scott boyd G stanford.edu). 


Materials Availability 


This study did not generate new unique reagents. 


Data and Code Availability 

Raw data from all serology Figures have been deposited on Mendeley at 
http://dx.doi.org/10.17632/j8r94pfrj6.1 and are publicly available as of the date of publication. АП 
original code has been deposited on the Zenodo platform at 
https://doi.org/10.5281/zenodo.5854880 and is publicly available at the Github repository 
https://github.com/boyd-lab/covid-infection-vs-vaccination as of the date of publication. Any 
additional information required to reanalyze the data reported in this paper is available from the 


lead contact upon request. 


EXPERIMENTAL MODELS AND SUBJECT DETAILS 

Plasma and saliva samples from Stanford BNT162b2 vaccinees 

To study antibody responses after first, second, and third dose vaccination with BNT162b2, we 
collected longitudinal blood and saliva samples from 59 vaccinees (29 were women, 27 were men, 


3 were unknown, and all donors were adults between the ages of 19 to 79 years). Baseline blood 
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samples were collected on day 0 before or immediately after the first vaccine dose. Individuals 
received their second dose on day 21 and a third dose about 9 months after the prime. Blood sample 
collection after prime vaccination was scheduled for days 1, 7, 21, 22, 28, and 42 with blood draws 
t one day from the assigned time point, or days 90, 120, or 210 € 1 week from the assigned time 
point. In addition, blood samples were collected оп days 0 to 3, 7 to 10, 21, and 28 after the third 
vaccine dose. Saliva samples were collected on day 42 after the first dose, as well as before, and 1 
or 2, 3 or 4, 7, and 21 days after the third vaccine dose. Peripheral blood was collected in vacutainer 
cell preparation tubes (CPTs) containing sodium citrate. Plasma was isolated and stored at -80°С. 
Saliva was collected from study participants, centrifuged, and supernatants were stored at -80?C. 
АП BNT162b2 vaccine study participants provided informed consent under Stanford University 


Institutional Review Board approved protocol IRB-55689. 


Plasma and saliva samples from Stanford COVID-19 patients 

Blood and saliva samples were collected between March and December 2020 from COVID-19 
patients who reported to Stanford Healthcare-associated clinical sites. SARS-CoV-2 infection was 
confirmed for all patients by RT-qPCR of nasopharyngeal swabs as described (Corman et al., 
2020; Hogan et al., 2020). Blood samples were collected in heparin- or EDTA-coated vacutainers. 
After centrifugation for collection of plasma or saliva, samples were stored at -80°C. The use of 
these samples for antibody testing was approved by the Stanford University Institutional Review 
Board (Protocols IRB-48973 and IRB-55689). 

Stanford COVID-19 patient cohort 1 included 530 plasma samples collected from 100 
moderately to severely ill COVID-19 patients (52 were women and 48 were men; ages ranged 


from 1 to 95 years; 23 were outpatients, 33 were admitted to hospital without needing ICU care, 
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20 were treated іп the ICU and 24 died of COVID-19) between March 2020 and August 2020. 
ELISA serology data for Wuhan-Hu-1 SARS-CoV -2-specific antibodies in these specimens have 
been reported previously (Róltgen et al., 2020). 

Stanford COVID-19 patient cohort 2 was included as a validation cohort of 87 samples from 74 
mostly mildly ill patients who had blood sample draws between March and December 2020, at 
approximately 21 days (n = 15 samples), 1 month (n = 23 samples), 3 months (n = 27 samples) 
and 7 months (n = 22 samples) after positive RT-qPCR testing for SARS-CoV-2 infection. Of 
those patients, 37 were women, 34 were men and 3 were unknown. Donors were 19 to 72 years of 
age and in terms of disease severity, 59 were mildly ill, 6 were moderate ill and 9 had a 
severe/critical disease course. Specimen time points were selected to match those of Stanford 
BNT162b2 vaccinee sample collections. Saliva samples were collected from five COVID-19 
patients. 

Stanford SARS-CoV-2 variant infection cohort blood samples were collected from COVID-19 
patients during acute infection with SARS-CoV-2 Alpha (n = 7) or Delta (n = 34) variants. Samples 
were from 20 women and 21 men, all between 2 and 92 years of age. SARS-CoV-2 genotyping 
data were obtained using a multiplex, mutation-specific RT-qPCR targeting N501Y, E484K, and 
L452R, as previously described (Wang et al., 2021). Samples from the first multiplexed reaction 
suspected to contain the Alpha variant were analyzed with a second confirmatory genotyping RT- 
qPCR assay to detect mutations encoding the N501Y amino acid change, as described (Dashdorj 


et al., 2021а). 
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Plasma samples from Mongolian vaccinees 

To study SARS-CoV-2 variant-specific IgG responses elicited by different COVID-19 vaccines, 
we tested plasma samples collected in July 2021 from 196 Mongolian vaccine study participants 
(109 were women, 87 were men, all were adults between 20 and 85 years of age) who had been 
fully vaccinated with one of four COVID-19 vaccines: BioNTech-Pfizer BNT162b2 (n = 47), 
AstraZeneca ChAdOx1-S (n = 50), Sputnik V Gam-COVID-Vac (n = 45) and Sinopharm BBIBP- 
CorV (n = 54). Participants were recruited by public announcement and volunteers were enrolled 
after signing the consent form approved by the Ethics Review Board at the Ministry of Health of 
Mongolia. SARS-CoV-2 pseudotyped virus neutralization and RBD-ACE2 blocking data on the 
same samples have been reported previously (Dashdorj et al., 2021a). Peripheral blood was 


collected in CPT, centrifuged for collection of plasma, and stored at -80°C. 


Healthy human control (HHC) plasma and saliva 

37 plasma and 20 saliva samples from HHCs collected before the onset of the COVID-19 
pandemic for studies at the Sean N. Parker Center for Allergy & Asthma Research were used to 
verify pre-pandemic antibody binding concentrations to the different coronavirus antigens, and 
manufacturer-provided cutoffs for positive serology assay results. Use of these samples was 
approved by the Stanford University Institutional Review Board (Protocols IRB-8629 and IRB- 


60171). No demographic information was available for these samples. 


Axillary LN core biopsies and post-mortem peribronchial LN tissues 


To analyze and compare GC architecture in response to COVID-19 vaccination and SARS-CoV- 


2 infection, we collected axillary LN core needle biopsies from BNT162b2 or mRNA-1273 
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vaccinees, and excised post-mortem peribronchial LNs from patients who died of COVID-19. For 
the selection of vaccinee tissues, we performed a retrospective search of our pathology archives 
and medical records between January 2021 and June 2021 for female patients who received either 
mRNA-1273 or BNT162b2 vaccination and subsequently underwent an ipsilateral axillary LN 
core needle biopsy due to mammographic findings and routine clinical care. Seven patients 
underwent biopsy one to eight weeks after vaccination with their second dose of mRNA vaccine. 
Three unvaccinated females undergoing axillary LN core biopsy for routine clinical care and 
mammographic findings served as controls. We included six peribronchial LNs from two female 
and four male patients who died of COVID-19 before August 2020, one to three weeks after 
symptom onset. Control post-mortem peribronchial LN biopsies were from pre-pandemic patients 
who died of non-COVID-19 causes. Autopsies were done by the Arizona Study of Aging and 
Neurodegenerative Disorders Brain and Body Donation Program (Beach et al., 2015). Analysis of 
these tissues was approved by Stanford University Institutional Review Board Protocol IRB- 


48973. 


METHOD DETAILS 

MSD ECL binding assays 

Plasma samples from vaccinees and COVID-19 patients were heat-inactivated at 56°С for 30 
minutes and tested using multiplexed ECL detection in a 96-well plate format with MSD? V- 
PLEX? serology panels and instrumentation according to the manufacturer's instructions. V- 
PLEX COVID-19 Coronavirus Panel 2 kits were used to detect IgM, IgG, and IgA antibodies to 
SARS-CoV-2 №, 51 NTD, КВР, and spike antigens and to spike proteins of SARS-CoV and other 


HCoVs including HCoV-OC43, HCoV-HKUI, HCoV-NL63, and HCoV-229E. V-PLEX SARS- 
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CoV-2 Panel 9 and 11 kits were used (о determine IgG antibody concentrations and RBD-ACE2 
blocking antibody percentages to different SARS-CoV-2 variant RBDs, with Alpha, Beta, Gamma, 
Epsilon, Kappa, Eta/Iota, B.1.526.2, P.3 and Wuhan-Hu-1 present in both panel 9 and 11, and with 
B.1.214.2 in panel 9 and Delta in panel 11. V-PLEX SARS-CoV-2 Panel 20 kits were used to 
determine IgG antibody concentrations to Alpha, Beta, Gamma, Delta, and Wuhan-Hu-1 SARS- 
CoV-2 variant spike proteins. Plasma samples were analyzed in duplicate at a 1:5,000 (for IgG 
binding assays) or a 1:10 (for RBD-ACE2 blocking assays) dilution in MSD diluent. Coronavirus- 
specific antibodies were detected with anti-human IgM, IgG, or IgA antibodies, or indirectly with 
human ACE2 protein (for RBD-ACE2 blocking assays) conjugated to SULFO-TAG™ ECL labels 
and read with a MESO? QuickPlex? SQ 120 instrument. Cutoff values for positive antibody test 
results for Wuhan-Hu-1 antigens were determined by the manufacturer based on sera from 200 
pre-pandemic healthy adults and 214 PCR-confirmed COVID-19 patients. We tested an additional 
37 healthy adult pre-pandemic plasma specimens to evaluate the manufacturer's cutoff values, and 
to determine cutoffs for positive binding to variant virus antigens, defined as the mean plus three 
standard deviations of the results from the pre-pandemic specimens. Antibody binding ratios for 
Wuhan-Hu-1 and viral variant antigens were only calculated for specimens that were above the 
cutoff values for positive results. Saliva samples were analyzed in duplicate at a 1:5 dilution in 
MSD diluent 2. Each plate contained duplicates of a 7-point calibration curve with serial dilution 
of a reference standard, a blank well and three positive control samples. Calibration curves were 
used to calculate antibody unit concentrations by backfitting ECL signals measured for each 


sample to the curve. 


MSD ECL spike antigen detection 
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SARS-CoV-2 spike antigen was quantified in plasma samples using an antigen capture ЕСІ, 
immunoassay platform (Meso Scale Discovery). S-PLEX? SARS-CoV-2 spike Kit assays were 
performed according to manufacturer instructions. A 7-point calibration curve and negative control 
consisting of assay diluent were run in duplicate on each plate. Plates were read using a MESO 
QuickPlex SQ 120 instrument. Raw signals were converted to a concentration based on linear 
regression to the 7-point calibration curve. Recombinant SARS-CoV-2 spike protein used for 


plasma spiking experiments was made by ATUM (https://www.atum.bio/). 


Histology, immunohistochemistry & in situ hybridization 

LN core needle and autopsy tissue samples were fixed in formalin and embedded in paraffin 
(FFPE), and sectioned. Once unstained slides were generated and initial Hematoxylin and Eosin 
(H&E) stained sections were analyzed, two distinct 0.6 mm areas from each LN sample were cored 
out of each tissue block and re-embedded to construct a tissue microarray (TMA). 
Immunohistochemistry was performed on four-micron sections using standard automated or 
manual methods including deparaffinization, peroxidase blocking, antigen retrieval, primary and 
secondary antibody incubation, detection with 3,3'-Diaminobenzidine (DAB) development, and 
counterstaining. Assays were performed on Roche Ventana (Tucson, AZ) Ultra instruments using 
Ventana Optiview detection, or Leica (Buffalo Grove, IL) Bond III instruments using Leica 
Polymer Refine detection or manually using Dako (Carpenteria, CA) Target Retrieval (TR) and 
Liquid DAB- Substrate Chromogen System with ImmPress (Vector, Burlingame, CA) secondary 
antibodies. 

For in situ hybridization, manual methods were used as previously described (Cloutier et al., 2021), 


using manufacturer-recommended protocols with the RNAScope 2.5 HD Assay-RED kit and 
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probes from Advanced Cell Diagnostics (Newark, СА). Two SARS-CoV-2 vaccine probes were 
developed to target bases 101-1143 of the spike encoding sequence of the BNT162b2 vaccine or 
bases 101-1488 of the spike encoding sequence of the mRNA-1273 vaccine. Both probes 
recognized SARS-CoV-2 mRNA vaccine, thus only the SARS-CoV-2 vaccine probe recognizing 
bases 101-1488 of mRNA-1273 vaccine are presented. To assess the specificity of SARS-CoV-2 
RNAScope vaccine probes, they were tested against SARS-CoV-2 infected placental tissue, in 
addition to staining for SARS-CoV-2 viral probe which targets bases 21631-23303 of the S-gene. 
SARS-CoV-2 RNAScope vaccine probes did not recognize SARS-CoV-2 virus. In addition, for 
each tissue tested there were internal negative control areas which did not react with SARS-CoV- 


2 vaccine probe. 


Whole Slide Imaging (WSI) and Quantitative Image Analysis 

Whole slide images of immunohistochemistry-stained slides for BCL6, PD-1 and CD21 were 
scanned at 40X magnification (0.25 uM per pixel) on an Aperio AT2 scanner (Leica Biosystems, 
Nussloch, Germany) in ScanScope Virtual Slide (SVS) format. Total lymphoid tissue and 
individual GCs were annotated as regions of interests (ROIs) using the QuPath open-source WSI 
software (Bankhead et al., 2017). GCs were defined as B cell areas with CD21- follicular dendritic 
cell networks and BCL6+ nuclei. Primary follicles were defined as B cell areas with CD21+ 
follicular dendritic cell networks without BCL6+ nuclei. For each GC ROI, the standard positive 
cell detection function was used to identify positive and negative cells with a single threshold of 
0.2 when scoring the cell DAB OD mean. For autopsy BCL6 analysis (COVID-19 and control) 
the threshold was adjusted to 0.05 given the dim expression of BCL6. For each ROI, the area in 


п>, number positive cells per mm”, and percent positive cells were calculated by QuPath. 
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A TMA slide containing two distinct 0.6 mm cores of each of the mRNA vaccine (n = 7) and 
vaccine control (n - 3) biopsies hybridized with SARS-CoV-2 mRNA RNAScope vaccine probe 
was scanned at 40X magnification (0.25 uM per pixel) on an Aperio AT2 scanner in SVS format. 
Total lymphoid tissue and individual GCs were annotated as ROIs using QuPath. For each GC 
ROI, the number of spots/clusters of RNAScope probe were detected using the QuPath subcellular 
detection option per manufacture instructions. Our detection parameters were as follows: 
Detection threshold - 0.6; split by intensity. Our split and cluster parameters were as follows: 


expected spot size = 2 um?, min spot size = 2 um”, max spot size = 3 um’. 


Co-detection by indexing (CODEX) 

All antibodies used for CODEX were first screened, titrated, and validated by individual staining 
on FFPE human tonsil tissue samples. Standard manual immunohistochemistry was used to cross- 
validate antibodies with the same, non-conjugated antibody clones listed in the key resource table. 
Tissue preparation for CODEX was undertaken by obtaining eight-micron thick sections from 
FFPE tissue blocks which were immobilized on charged square glass coverslips coated with 
polylysine (Electron Microscopy Sciences, Hatfield, PA) prepared according to the manufacturer's 
instructions. The coated glass coverslips were stained with a cocktail containing nucleotide- 
barcoded primary antibodies. The coverslips underwent nuclear staining (DAPI) and were loaded 
on the stage of an automated inverted fluorescence microscope connected to the robotic fluidic 
system known as CODEX (Akoya biosciences, Marlborough, MA). In addition to the nuclear stain, 
fluorophore-tagged complementary nucleotide sequences (reporters) were used to iteratively 
reveal three antibodies at a time per cycle. Two additional blank cycles at the beginning and at the 


end of the antibody reveal cycles were added for purposes of subtraction of auto-fluorescence 
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background. Automated image acquisition, processing, segmentation, апа fluidics exchange were 
conducted using an Akoya CODEX instrument and CODEX driver software (Akoya biosciences) 
(Black et al., 2021; Goltsev et al., 2018; Schürch et al., 2020). Composite images providing 
multiplex fluorescent signals of stained cell types, cellular niches, and tissue architecture were 


captured using CODEX® МАУ and FIJI software (Akoya biosciences). 


QUANTIFICATION AND STATISTICAL ANALYSIS 

Statistical tests were performed in R using base packages for statistical analysis and the ggplot2 
package for graphics. Box-whisker plots show median (horizontal line), interquartile range (box), 
and the end of the lower whisker representing the smallest observation greater than or equal to the 
2596 quantile minus 1.5 times the interquartile range, and the end of the upper whisker representing 
the largest observation less than or equal to the 75% quantile plus 1.5 times the interquartile range. 
In serological analyses where statistical significance was tested, significance was defined as: ***p 
value < 0.001; **p value < 0.01; *p value < 0.05. 

For the principal component analysis (PCA) we log-transformed, calculated z-scores, and ran PCA 
on MSD antibody concentration measurements or Wuhan-Hu-1/variant RBD IgG concentration 
ratios from a reference time point after COVID-19 vaccination or SARS-CoV-2 infection using 
Python v3.7.10 and packages numpy v1.19.1, pandas v1.2.5, and scikit-learn v1.0. We then applied 
these transformations to matching data from all other time points, enabling us to visualize the 
change over time in these serology measurements on a consistent PCA reference. Plots were 


created with Python packages matplotlib version 3.3.2 and seaborn version 0.11.2. 
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For the analysis of the homogeneity or dispersion of serology measurements in groups differing 
by vaccination or infection status (Figures 2C and 2D) at a particular time point, we plotted each 
group's distribution of Euclidean distances to its centroid (calculated with Python package scipy 
version 1.6.2). These distance distributions were consistent when calculated in the raw 
measurement data space of arbitrary units (AU) for the MSD ECL assay or in the transformed PC1 
and PC2 space after embedding the raw measurements into the PCA space created from one 


reference time point. 


Ratios of concentrations of IgG binding to Wuhan-Hu-1 RBD compared to variant virus RBDs 
were plotted for specimens with IgG binding above the cutoff for positive binding to Wuhan-Hu- 
1 RBD (Figures 3, 4 and 54) to avoid distortion of ratios by samples without specific binding. 
Ratios of IgG binding to spike antigens were calculated in a similar manner (Figure SAC, lower 
panel). Corresponding IgG concentrations for samples used to calculate ratios were plotted for 
reference (Figures ЗА, 4А, 54А and S4C upper panel). To quantify serological imprinting from 
prior Wuhan-Hu-1 antigen exposure from vaccination on subsequent responses to breakthrough 
infection with the Delta variant, we first computed the ratio of Wuhan-Hu-1 RBD binding level to 
Delta RBD binding level. Here, a ratio of one indicates even preference, while ratios greater than 
one indicate preferential binding of Wuhan-Hu-1 over Delta. Each ratio is symmetric with its 
inverse; for example, а Wuhan-Hu-1/Delta binding ratio of 4/5 indicates the same degree of 
preference for Delta binding as the ratio 5/4 indicates for Wuhan-Hu-1 binding preference. We 
then log-transformed the ratios, which sets the even preference level at zero, with positive values 
corresponding to Wuhan-Hu-1 preference, and makes these values symmetric around zero (e.g., a 


value of -0.2 indicates the same level of preference for Delta as 0.2 does for Wuhan-Hu-1 


41 


921 


922 


923 


924 


925 


926 


927 


928 


929 


930 


931 


932 


933 


934 


935 


936 


937 


938 


939 


940 


941 


942 


943 


binding). Finally, we rescaled the negative and positive values separately to the ranges -100% to 
0 and 0 to +100%. The resulting magnitudes are binding preferences relative to the maximum 
binding preference observed for a particular variant, including data from individuals only exposed 
to Wuhan-Hu-1 or individual variants. In particular, -100% refers to the maximum observed 
binding preference towards Delta. We plotted the distributions of these binding preference levels 
for BNT162b2 vaccinees without known prior SARS-CoV 2 infection, for Delta infection cases 
with no recorded prior infection or vaccine exposure, and for Delta breakthrough infections 
following Wuhan-Hu-1 wild-type vaccination. An example binding preference level of +30% for 
a Delta breakthrough infection case suggests this individual is at 3096 of the most imprinted state 
for Wuhan-Hu-1 preference. А binding preference of -20% would instead suggest that, following 
infection, this individual lost imprinting to Wuhan-Hu-1 and gained preference for binding the 


Delta variant. 


Supplemental information titles and legends 

Figure S1. Anti-SARS-CoV-2 Ig antibody responses in plasma and saliva following 
BNT162b2 vaccination, related to Figure 1. 

Anti-SARS-CoV-2 N, RBD, and spike (S) IgG (A), IgM (B), and IgA (C) responses are shown for 
plasma from individuals who received BNT162b2 prime (DO, n = 59) and second dose (D21, n = 
58) vaccination. Box-whisker plots of the WHO binding arbitrary unit (BAU/mL) anti-SARS- 
CoV-2 concentrations show the interquartile range as the box and the whisker ends as the most 
extreme values within 1.5 times the interquartile range below the 2596 quantile and above the 7596 
quantile. Comparisons between groups of previously SARS-CoV -2-infected (CoV-2+) versus non- 


infected individuals and female versus male were by the two-sided Wilcoxon rank sum test; 
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comparison between age groups (< 40; 40 to 60; > 60 years) was done using pairwise Wilcoxon 
rank sum test with Bonferroni correction. *p « 0.05, **p « 0.01, ***p « 0.001. 

(D) Anti-SARS-CoV-2 М, RBD, and 5 IgG concentrations in BAU/mL are shown for saliva from 
individuals who received BNT162b2 prime/boost and 3 dose vaccination (upper left panel). Anti- 
SARS-CoV-2 М, RBD, and 5 (upper right panel) concentrations in BAU/mL, as well as anti- 
SARS-CoV-1 and anti-HCoV-OC43, -НКОІ, -NL63, and -229E S IgG (lower panel) 
concentrations in MSD arbitrary units (AU/mL) are shown for saliva collected on D42 after 
BNT162b2 prime vaccination (vaccinee), around D42 post-symptom onset for COVID-19 patients 
(CoV-2+), and before the onset of the COVID-19 pandemic for pre-pandemic healthy human 
controls (Pre-pan). Box-whisker plots of anti-SARS-CoV-2 IgG concentrations show the 
interquartile range as the box and the whisker ends as the most extreme values within 1.5 times 
the interquartile range below the 25% quantile and above the 75% quantile. Statistical test for 
significance between groups (СоУ-2-; Pre-pan, Vaccinee) was performed using pairwise 


Wilcoxon rank sum test with Bonferroni correction. *p < 0.05, **p « 0.01, ***p « 0.001. 


Figure S2. The magnitude of antibody responses is not correlated with reported vaccine- 
associated side effects (SEs), related to Figure 1. 

(A) Frequency of site-specific and systemic vaccine-associated SEs after prime (light green) and 
second dose (dark green) BNT162b2 vaccination. 

(B) Box-whisker plots of the MSD AU/mL anti-SARS-CoV-2 IgG concentrations in BNT162b2 
vaccinee plasma collected оп D28 post-vaccination show the interquartile range as the box and the 
whisker ends as the most extreme values within 1.5 times the interquartile range below the 2596 


quantile and above the 7596 quantile. For a given SE (rows), vaccinees were grouped according to 
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SEs were unknown are shown as white boxplots. 


Figure S3: BNT162b2 vaccination produces less broad serological responses to endemic 
HCoVs compared to SARS-CoV-2 infection, related to Figure 2. 

(A, B) Anti-SARS-CoV-1 spike, and anti-HCoV-OC43, -HKUI, -NL63, and -229E spike IgM, 
IgG, and IgA antibody responses are shown for individuals who received BNT162b2 prime (DO) 
and boost (D21) vaccination doses and for COVID-19 patients. 

(A) The heatmap shows the development of antibody responses in longitudinal samples from 
vaccinees/patients collected at / during DO, D7 / week 1, D21 / weeks 2&3, D28 / week 4, D42 / 
weeks 5&6, and D90 / >week 7 after vaccination / COVID-19 symptom onset (x-axis). The color 
scale encodes the median values of log10 MSD AU/mL concentrations. 

(B) Box-whisker plots show the development of antibody responses in longitudinal samples from 
vaccinees / patients collected at / during DO, D7 / week 1, D21 / weeks 2&3, 028 / week 4, D42 / 
weeks 5&6, and D90 / >week 7 after vaccination / COVID-19 symptom onset (x-axis). Box- 
whisker plots show the interquartile range as the box and the whisker ends as the most extreme 
values within 1.5 times the interquartile range below the 25% quantile and above the 75% quantile. 
Statistical test: pairwise Wilcoxon rank sum test with Bonferroni correction. *p < 0.05, **p « 0.01, 
***p < 0.001. 

Individuals were classified as vaccinees who have not been previously exposed to SARS-CoV-2 
(Vaccinees); outpatients (Outpt) and hospital admitted patients (Admit); intensive care unit (ICU) 


patients, and those who died from their illness (Death). 
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(С) Box-whisker plots show anti-SARS-CoV-1 spike and anti-HCoV spike antibody responses in 
plasma samples from individuals who received BNT162b2 prime (DO, n = 59 individuals), second 
dose (D21, n -58 individuals) and third dose (around month 9, n - 36 individuals) vaccination. 
Box-whisker plots show the interquartile range as the box and the whisker ends as the most extreme 


values within 1.5 times the interquartile range below the 25% quantile and above the 75% quantile. 


Figure S4: Greater breadth of IgG binding to SARS-CoV-2 variant RBDs following 
BNT162b2 vaccination compared to infection with Wuhan-Hu-1 SARS-CoV-2 (Validation 
cohort), related to Figure 3. 

(A, B) Anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant RBD IgG responses are shown for 
Stanford individuals who received BNT162b2 vaccination and for Wuhan-Hu-1-infected COVID- 
19 Stanford patient cohort 2 at different time points after vaccination / COVID-19 symptom onset. 
Box-whisker plots show the interquartile range as the box and the whisker ends as the most extreme 
values within 1.5 times the interquartile range below the 25% quantile and above the 75% quantile. 
Significance between groups were tested with pairwise Wilcoxon rank sum test with Bonferroni 
correction. *p < 0.05, **p « 0.01, ***p « 0.001. 

(A) Anti-RBD IgG concentrations. 

(B) Ratios of anti-Wuhan-Hu-1 to variant RBD IgG concentration. 

(C) Anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant spike IgG responses as anti-spike IgG 
concentrations (upper panels) and as ratios of anti- Wuhan-Hu-1 to variant spike IgG concentration 
(lower panels) are shown for Stanford individuals who received BNT162b2 vaccination and for 
Wuhan-Hu-1-infected COVID-19 Stanford patient cohort 1 and 2 samples. Box-whisker plots 


show the interquartile range as the box and the whisker ends as the most extreme values within 1.5 


45 


1012 


1013 


1014 


1015 


1016 


1017 


1018 


1019 


1020 


1021 


1022 


1023 


1024 


1025 


1026 


1027 


1028 


1029 


1030 


1031 


1032 


1033 


1034 


times the interquartile range below the 25% quantile and above the 75% quantile. Significance 
between groups were tested with pairwise Wilcoxon rank sum test with Bonferroni correction. 
*p < 0.05, **p « 0.01, ***p « 0.001. 

(D) Percentage blocking of ACE2 binding to RBD of specified viral variants by plasma 


antibodies of BNT162b2 vaccinees and Stanford patient cohort 2 samples. 


Figure S5: Anti-SARS-CoV-2 RBD IgG signatures following BNT162b2 vaccination and 
SARS-CoV 2 infection, related to Figure 5. 

(A) Ratios of anti-Wuhan-Hu-1 to variant RBD IgG concentration are shown for Stanford 
individuals who received BNT162b2 vaccination at different time points after second dose (D21, 
n =58 individuals) and third dose (around month 9, n = 36 individuals) vaccination. Box-whisker 
plots show the interquartile range as the box and the whisker ends as the most extreme values 
within 1.5 times the interquartile range below the 25% quantile and above the 75% quantile. 

(B) Principal component analysis (PCA) of anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant RBD 
IgG concentrations across Stanford BNT162b2 vaccinees, Stanford COVID-19 patient cohort 2 


and SARS-CoV-2 variant-infected patients. 


Figure S6: Disrupted LN GCs in COVID-19 patients versus mRNA vaccinees, related to 
Figure 6. 

(A) LN GC histology for COVID-19 patients (left) and mRNA vaccinees (right) evaluated with 4- 
color Co-detection by indexing (CODEX) immunofluorescence analysis for CD20 (red), CD3 
(blue), BCL6 (magenta) and CD21 (yellow) markers of B cells, T cells, GC B cells (or T follicular 


helper cells) and follicular dendritic cells, respectively. 
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(B) Representative CD21 immunohistochemistry of secondary (left) and primary (right) follicles 


of four autopsy patients who died of COVID-19 and two control autopsy patients. 
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Highlights 


е Vaccination confers broader IgG binding of variant RBDs than SARS-CoV-2 
infection 

e Imprinting from initial antigen exposures alters IgG responses to viral variants 

е Histology of mRNA vaccinee lymph nodes shows abundant germinal centers 

e Vaccine spike antigen and mRNA persist for weeks in lymph node germinal 
centers 


In Brief 

Human antibody responses to SARS-CoV-2 differ between vaccination and infection, 
with vaccination (regardless of vaccine type) inducing more productive lymph node 
germinal center responses and a broader range of IgG neutralizing antibodies capable 
of recognizing viral variants. 
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